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EFFECT  OF  NONUNIFORM  SIZE  ON  INTERNAL 
STRESSES  IN  A  RAPID,  SIMPLE  SHEAR  FLOW  OF 
GRANULAR  MATERIALS 
Part  2.  Multiple  Grain  Sizes 

Hayley  H.  Shen 


INTRODUCTION 

The  flow  of  a  granular  material  is  a  dynamic  phenomenon  that  can  be  modeled  by  macroscopic 
equations  describing  the  balance  of  mass,  momentum  and  energy  within  the  system.  A  correct 
formulation  of  the  stresses  and  the  rate  of  energy  input-output  is  crucial  to  the  success  of  a  math¬ 
ematical  model  of  such  a  system.  The  mathematical  formulations  of  the  stresses  and  the  rate  of 
energy  input-output  are  called  the  constitutive  equations.  The  constitutive  equations  for  a  flow¬ 
ing  granular  material  can  range  over  a  wide  spectrum,  from  Newtonian  to  highly  non-Newtonian 
forms. 

Although  a  nondimensional  parameter  that  determines  the  form  of  the  constitutive  relations 
is  unknown  yet,  it  is  believed  that  the  flow  rate,  the  concentration  and  the  material  properties  of 
the  granular  solids  determine  the  regime  of  the  granular  flow  and  hence  the  form  of  the  constitu¬ 
tive  equations. 

Within  a  given  flow  regime,  the  constitutive  relations  are  determined  by  the  way  energy  and 
momentum  are  transferred  microscopically  from  one  location  to  another.  Different  transferring 
mechanisms  produce  different  constitutive  equations.  Two  examples  of  such  different  constitu¬ 
tive  equations  can  be  found  in  Bagnold  (1954)  where  either  a  linear  or  a  quadratic  stress-strain- 
rate  relation  is  used,  depending  on  the  flow  regime,  in  slow  flow  viscous  dissipation  in  the  granu¬ 
lar  interstice  dominates  and  gives  a  linear  stress-strain-rate  relation,  while  in  a  fast  flow,  dissipa¬ 
tion  by  collisions  among  granules  dominates  and  gives  a  quadratic  stress-strain-rate  relation  simi¬ 
lar  to  those  in  a  turbulent  fluid  flow.  Another  type  of  constitutive  relation  has  been  derived  for 
deforming  ice  fields  in  the  polar  seas  where  pressure  ridge  formation  was  the  major  dynamic  acti¬ 
vity  (e.g.  Rothrock  1975).  It  is  found  that  the  stress-strain-rate  relation  is  of  the  plastic  type. 
Recently,  however,  there  have  been  some  doubts  about  the  applicability  of  this  constitutive  law 
to  the  ice  edge  where  concentration  is  too  low  to  form  ridges. 

In  the  past  few  years,  a  number  of  researchers  have  formulated  constitutive  equations  for  granu¬ 
lar  flows  (Kanatani  1979,  Ogawa  et  al.  1980,  Shen  and  Ackermann  1982,  Campbell  and  Brennen 
1983,  Walton  1983,  Jenkins  and  Savage  1983  and  Lun  et  al.  1984),  and  most  of  these  have  as¬ 
sumed  intergranular  collisions  to  be  the  major  momentum  transfer  and  energy  dissipation  mechan¬ 
ism.  The  resulting  constitutive  equations  are  quadratic  and  check  well  with  experimental  data 
(Shen  and  Ackermann  1982).  The  above  researchers  found  that  the  stress  level  generated  in  such 


a  granular  tlow  depends  crucially  on  material  properties  of  the  physical  constituents:  the  restitu¬ 
tion  and  frictional  coefficient,  densities  of  the  granular  solid  and  the  interstitial  fluid,  and  the 
drag  coefficient  of  the  tluid. 

While  able  to  give  a  great  deal  of  physical  insight  into  the  dynamics  of  a  fast  granular  flow,  the 
above  studies  assumed  that  the  granular  solids  are  either  disks  or  spheres  and  are  uniform  in  size. 
Natural  materials  of  interest  are  generally  irregular  in  shape  and  have  various  sizes.  As  evidenced 
by  experiments  on  sediment  transport  (e.g.  Gilbert  1914,  Durand  1953),  under  the  same  driving 
mechanism,  a  widely  graded  sediment  is  transported  at  a  much  higher  rate  than  a  relatively  uniform¬ 
ly  graded  sediment.  These  phenomena  suggest  that  the  stress  levels  and  the  rate  of  energy  dissipa¬ 
tion  are  significantly  influenced  by  the  size  distribution  existing  in  the  flowing  granular  material. 

This  report  presents  the  formulation  of  constitutive  equations  incorporating  a  size  distribution 
in  the  granular  material.  Intergranular  collision  is  assumed  the  major  dynamic  activity  at  the  micro¬ 
scopic  level,  as  in  an  earlier  work  where  granular  materials  of  one  size  were  considered  (Shen  and 
Ackeimann  1982).  Because  of  the  present  limited  knowledge  of  how  to  treat  shape  effects,  the 
analysis  will  be  confined  to  the  How  of  either  disks  or  spheres.  Nevertheless,  the  result  of  this 
work  provides  information  necessary  to  a  more  realistic  analysis  of  virtually  any  natural  granular 
flow,  including  drifting  ice  fields  in  the  polar  seas  and  avalanches  of  snow  or  debris. 

Die  analysis  for  a  complete  spectrum  of  sizes  is  based  on  the  result  of  analyzing  a  two-size  mix¬ 
ture  ( Shen  1 985 ).  f  or  the  sake  of  completeness,  the  analysis  for  a  two-size  mixture  is  reviewed 
briefly  in  what  follows;  the  details  are  given  in  the  companion  to  this  report,  Effect  of  nonuniform 
size  mi  internal  stresses  in  a  rapid,  simple  shear  /low  of  granular  materials :  Part  1.  Two  grain  sizes. 

Hie  two-size  mixture  considered  here  is  an  assembly  of  spherical  particles  because  spherical 
particles  model  a  wide  variety  of  natural  materials  and  because  the  only  existing  experimental 
data  for  direct  comparison  with  the  theory  were  obtained  for  spherical  particles.  However,  it  is 
a  st raiglit forward  extension  to  apply  the  same  analysis  to  disks.  The  complete  spectrum  analysis 
will  be  carried  out  for  both  spheres  and  disks. 


CONSTITUTIVE  EQUATIONS  FOR  A  TWO-SIZE  MIXTURE 

Assume  a  homogeneous  mixture  ot 
spheres  with  two  sizes  under  a  simple  shear 
motion  as  illustrated  in  f  igure  1 .  The  in¬ 
terstitial  tluid  effect  is  not  included  in  this 
model.  In  a  crude  way.  the  interstitial  fluid 
effect  can  be  treated  as  did  Ackermann  and 
Shen  ( 1 982 ).  Die  diameters  of  the  large 
and  small  spheres  are  Dv  and  Ds  respective¬ 
ly.  Although  extending  this  model  to  mix¬ 
tures  of  different  material  properties  is  a 
st  raiglit  forward  matter,  in  the  present  studs 
only  the  size  effect  is  considered.  The  ma¬ 
terial  property  tor  all  spheres  is  thus  the 
same  and  includes  the  density  ps.  the  resti-  Figure  I.  Shear  flow  of  a  mixture  of  spheres. 

tution  coe'ficient  e.  a  rid  the  frictional  co¬ 
efficient  p. 

l  et  <7,  and  <7S  denote  the  number  percentages  of  the  large  and  small  spheres  in  a  unit  volume 
respectivek  .  This  means  that  it  <7,  :<7's  =  7()  30.  there  are  70  large  spheres  and  30  small  ones  in 
a  random  sample  ot  100  sphetes.  I  he  ratio  ot  large  to  small  spheres  surrounding  an  average  sphere 
is  1*  |  r,s.  I  or  a  given  combination  ot  (7j  <7S  and  f\  .  let  C  n  denote  the  volume  concentra¬ 
tion  ot  the  spheres  when  all  neighboring  spheres  aie  touching,  i.e..  C  n  is  the  volume  concentration 
tor  the  densest  random  packing.  It  all  spheres  are  then  separated  by  the  distance  s  away  from  the 
neighboring  ones,  the  volume  concentration  (  tor  this  looser  state  is 


(1) 


C=C0[R/(R+s)}3 

where  R  is  the  radius  of  a  typical  “cell”  formed  by  a  given  sphere  and  all  its  neighboring  spheres 
(the  power  becomes  2  for  disk -shaped  material),  and  is  given  as  (Shen  1985): 

r  =  glr1  +  csrs  (2) 

where 


RL  -  ('l.Dl 

+  Gs</)l  +  Ds  )/2 

(3a) 

Rs  =  Gl(D, 

+  n s  >/:  +  c;s/)s . 

(3b) 

When  sheared,  the  assembly  of  these  spheres  will  start  to  collide  with  each  other  and  produce 
velocity  fluctuations.  Hie  momentum  transfer  caused  by  these  collisions  gives  rise  to  stresses 
which  are  modeled  as  (Bagnold  1954): 

hi=/V  f  ■  ^  (4) 

where  Ttj  =  stress  in  the  Jtj  direction  on  a  surface  normal  to  the  .tj  direction 
A M  -  momentum  transferred  in  the  x^  direction 
f  =  collision  frequency/2 

p{  =  number  of  spheres  per  unit  area  normal  to  the  jo  direction. 

The  formulation  of  stresses  expressed  in  eq  4  does  not  incorporate  the  fact  that  there  are  two 
kinds  of  spheres,  in  which  case  we  consider  the  control  volume  shown  in  Figure  2.  Stresses  are 
produced  by  four  different  kinds  of  collisions.  The  correct  form  of  eq  4  for  a  two-size  mixture 
is  thus 


=  L 
P.Q 


AM, 


^PQ 


>’Qi 


P  Pi 


(5) 


e  A MPq.  is  the  momentum  transfer  in  the  direction  when  a /’-size  and  a  0-size  sphere 
collide  ( P.Q  can  be  large  or  small  as  denoted  by  L  or  S),  A^,q  is  the  frequency  of  collisions  any 
given  /"-size  sphere  receives  from  0-size  particles,  and  ppj  is  the  number  of  P- size  spheres  on  a 
unit  area  normal  to  the.tj  direction. 


Q  ,9  Exterior  Porticles 
Surfoce  Porticles 
0,0  Interior  Particles 


Figure  2.  A  control  volume  in  a  granular  flow  of  uniform  spheres. 


shear  stress  (two-si/e  with  spectrum  given  by  Figure  3b)  _  1 

shear  stress  (single-si/e  with  spectrum  given  by  Figure  3a)  2.5 

It  is  worthwhile  noting  that  in  all  the  stress  computations,  Cu=  0.74  was  used.  This  value  of  CD 
corresponds  to  the  hexagonal  dense  packing  of  single-size  spheres.  In  a  granular  material,  the  dens¬ 
est  state  in  which  all  neighboring  spheres  are  touching  may  not  exhibit  the  hexagonal  structure. 
Moreover,  the  existence  of  more  than  one  size  in  the  flowing  assembly  can  further  vary  the  value 
of  C\t.  The  effect  of  size  distribution  on  the  value  of  C0  has  been  studied  (e.g.  McGaw  1967, 
Visscher  and  Bolsterli  1972),  but  is  not  sufficiently  understood  for  the  purpose  of  this  work.  Fur¬ 
ther  study  is  needed  to  incorporate  the  variability  of  Cu  as  a  function  of  the  size  distribution. 

As  an  example,  the  results  shown  in  eq  50a  and  56a  are  applied  to  the  log-normal  distributions 


/'(/>)=  :  exp  -  -  (  —  )  I 

x/2/r  oD  '  2  a  / 


(57) 


The  result  is  shown  in  Figure  8  with  77  =  0, 
n min  =  D nlax  =  er’  +  3°  and 


f  =  r2 !  (log-normai )/r2 1  ( />M  ) 


where  Ds.  is  the  mean  size  of  the  log-normal 
distribution  and  the  shear  stress  r21  (Z)M )  is 
computed  assuming  a  single-size  assembly 
(Shell  and  Ackermann  1982,  1984). 


CONCLUSION 


Theoretical  analysis  of  the  stress-strain-rate  Figure  8.  The  effect  of  log-normal  size  dis- 

relationships  for  rapid  granular  flows  have  been  trihution  on  the  shear  stresses  in  a  simple 

limited  to  single-size  assemblies  of  spheres  or  shear  flow  o)  spheres  or  disks. 

disks.  Fxtension  to  an  assembly  with  an  arbi¬ 
trary  size  distribution  is  needed  for  analyzing 

real  granular  materials.  A  study  was  carried  out  for  a  two-size  mixture  of  spherical  particles  (Shen 
1985 ).  In  that  study ,  assumptions  were  made  that  restrict  the  applicability  of  the  result  to  cases 
when  the  size  ratio  is  nearly  1 . 

Two  studies  are  made  in  this  report.  First,  the  stress-strain-rate  relationship  is  established  for 
a  binary  spherical  assembly  at  the  other  extreme,  namely,  as  the  ratio  of  the  two  sizes  is  nearly  0. 
With  the  combination  of  the  two  extreme  cases  when  the  size  ratio  is  nearly  1  or  0,  a  smoothing 
technique  is  used  to  give  an  approximate  stress  state  for  any  intermediate  size  ratio.  Flence,  the 
constitutive  equations  for  a  binary  assembly  of  spheres  of  an  arbitrary  size  ratio  are  obtained. 
Second,  an  arbitrary  spectrum  of  size  distribution  is  incorporated  into  the  analysis.  The  effect  of 
size  distribution  of  the  stresses  is  quantified  for  both  the  sphere-shaped  and  the  disk -shaped  parti¬ 
cles.  The  result  demonstrates  a  reduction  of  stresses  as  the  size  gradation  broadens.  This  phenom¬ 
enon  may  explain  some  of  the  observations  from  earlier  experimental  work  on  sediment  transport. 

As  indicated  earlier  in  this  report,  further  work  is  needed  to  improve  the  derivation  for  the  av¬ 
erage  gap  between  neigh  boring  spheres,  since  the  densest  solid  concentration,  although  assumed 
here  to  be  a  constant,  is  a  function  of  the  size  distribution.  This  is  true  for  both  mixtures  of  two 
sizes  and  of  a  complete  size  spectrum.  The  explicit  form  of  this  function  is  unknown  yet. 

Also,  in  the  case  of  a  complete  'occtruin  of  sizes,  as  the  size  gradation  widens,  the  larger  parti¬ 
cles'  fluctuations  reduce  to  zero  because  of  the  equipartition  of  energy.  In  this  case,  the  analysis 
has  to  he  modified  to  reevaluate  the  dynamic  interactions  involving  larger  particles.  For  a  binary 
mixture,  this  modification  has  been  carried  out  in  this  report.  However,  it  is  not  clear  how  to  pro¬ 
ceed  along  this  line  for  a  general  spectrum  of  sizes.  Nevertheless,  based  on  the  result  of  a  binary 


Psd ,,  ,  ,  du  /£>+£>'  ,  \  D2D'2 

6  (1+e)^2(  2  7  D2+D'2 

(52a) 

p  d  t -  n2  n'2 

^I)D'2  “  ^  0+e)\AD  +  VD'  £)2+q'2 

(52b) 

_  .  nd  /  1-e2  p(l+6)  P2(l+e)2\  ,  n  D2d'7 

-  Ps  4  U  +  n  4  )  (1D  '>D2W'2  ■ 

(52c) 

The  equipartition  of  energy  becomes 

D2v\t=D'2v'2.  . 

(53) 

Substituting  eq  41  d,  51a-d,  52a  and  c  into  eq  47,  we  obtain 

v2D  = 

[» 


<>’//(¥-) 


Cr>C|''  IrfOVD 


Z)2+D'2  £>' 


r rcDcu 
JJ  Z)3D'3 


where 


(54) 


A-  = 


(55) 


The  shear  and  normal  stresses  are  computed  by  substituting  eq  41  d,  51a  and  c  and  52a  and  b  into 
eq  39: 


T2l  =  Ps 


1+e  du 
3? r  t/jf 2 


ff/D+D'  \  D  „  „ 

JJ\  2  '  +S)  s  CDCo'  p 


^  dD’dD 
2+Z?  2  D 


r^Dvb 
J  D2 


dD 


(56a) 


*22  =  -Ps 


2(1+6) 


(56b) 


where  r>'()  is  defined  in  eq  54  and  s  is  defined  in  eq  5 Id. 

When  P(D)  equals  the  Dirac  delta  function  5(D),  eq  50a  and  b  and  56a  and  b  reduce  to  the 
stresses  for  a  single-size  sphere  or  disk  assembly  as  given  by  Shen  and  Ackermann  (1982*,  1984). 

If  eq  50a  is  used  for  a  complete  spectrum  of  sizes  to  reevaluate  the  comparison  between  the 
theory  and  the  experimental  data  obtained  by  Savage  and  Sayed  (1980)  discussed  previously,  the 
new  theory  gives 


•Because  of  an  underestimation  of  two  particles’  relative  fluctuation,  the  calculated  stresses  for  spheres 
should  have  been  l/y/2  of  what  were  shown  by  Shen  and  Ackermann  (1982). 
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If  eq  41-46  are  substituted  into  eq  47,  the  fluctuation  speed  for  the  mean  size  DM  is  obtained: 


>2  ..... 

‘'dm  “ 


a(fe>7/(er  * 


n 3 


ffc°c ■>• 


( DD’T */2  dD'tlD 


where 


(48) 


k  = 


(1+e)  (0,05  +  0.08  tx) 

1-c2  juO+c)_  u2(\+e)2 
4  n  4 


(49) 


For  any  given  size  between  D  and  D+dD,  the  fluctuation  speed  can  be  derived  using  eq  46  and  48. 

After  substituting  eq  41c  and  d,  43a  and  b,  46  and  48  into  eq  39,  the  shear  stress  r2i  and 
normal  stress  r22  in  a  simple  shear  flow  of  spheres  are 


r21  =  ps(  1  +e)  (0.05  +  0.08  p)  ~- 


rrn+n'  n  , 

ff(—  +s)jCDCD'  D3+D'3  dDt 


dD 


(50a) 


If? 


D  ( D/D'f 2  CDCD' 


Tji  =-PS^ 
7 r 


2<iW  *  DZ  D 

f^’  'o«D 


dD'dD 


(50b) 


These  were  computed  with  a  PC  Zenith-100  using  a  program  written  in  PASCAL.  The  computa¬ 
tion  only  takes  a  few  seconds  to  run  for  all  the  distribution  functions  tested. 

Based  on  the  analysis  of  a  planar,  simple  shear  flow  of  a  single-size  disk  assembly  (Shen  and 
Ackermann  1984),  the  above  analysis  for  a  spherical  assembly  can  easily  be  extended  to  a  disk¬ 
shaped  material  with  a  uniform  thickness.  Consider  a  layer  of  disks  with  uniform  thickness  d  and 
a  distribution  of  diameters  as  given  in  Figure  7a  or  7b.  This  layer  of  disks  is  undergoing  a  simple 
shear  flow  in  its  own  plane.  Equation  39  can  be  applied  directly  to  evaluate  the  stresses  generated 
by  disk -disk  collisions.  Using  the  notations  defined  before,  this  time  for  a  disk-shaped  assembly, 


CpID2 

0  fCuID2dD 


- 


4  Or 


D2d 


P 

[)  ir  Dd 


C„  V4 

s  =  /?<  f)  -R 


where  R  is  defined  in  eq  42b,  and 


(51a) 

(51b) 

(51c) 

(5  Id) 
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(41a) 


Gf,  - 


jV, 


U  fCtf/D3  dD 

6  Q) 


7T  Z)3 


=  - 


6  CD 


TT  DL 


i  t 

. ,s  (dr>vn 

DD'  =  fG^JD  '  s 


where 


C  1/3 

s  =  R(~)  -  R 


and 


R  =//gdG'd 


<»£  dD'dD  . 


The  average  momentum  transfer  and  energy  dissipation  are 

...  ,,4,  ?_  tr  .  .  .  du  ,D+D'  .2D3 D'3 

=  ^dd'  *  '  "  p*  “6  (1+e>  (°  05  +  0  08  d)  ^  (—  +  *) 

^Dl >'2  =  ^DD'  */  =  “f’s  3^  (1+e)V^D  +  ^D'  £)3  +£,'3 

,  _  ,  ZLfl-e2  .  P(l+e)  /a2(l+e)2  ,2  ,2  D3D'3 

*'»>’  -  ^  6  (  4  +  “V"  •  4  }  +  V)  W 

The  equipartition  of  energy  takes  the  same  form  as  before: 

DVd2  =  D'3Vp> . 


(41b) 

(41c) 

(4  Id) 


(42a) 

(42b) 

(43a) 

(43b) 

(43c) 

(44) 


Given  an  arbitrary  size  distribution  designated  by  some  probability  density  function  P(D),  the 


mean  size  0M  is 


Dm  =  /  DP(D)dD  . 

The  fluctuation  speed  for  this  mean  size  is 


'2  #2  /  Pd  3 

‘Tjm  =  l’t)  (^-3 


M 


(45) 


(46) 


The  energy  balance  equation,  eq  9,  now  takes  the  form 


du 
dx  2 


//' 


Ns  > 
/vnn 


f1'  P0dD'dD  =  JfEnn.  Nd  ^  dD’dD  ■ 


(47) 
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a.  Probability  distribution. 


b.  Probability  density. 


Figure  7.  Particle  size  distribution. 


N5  ' 

AMdd'  dD'  PDdD  .  (38) 

Hence  the  total  stress  caused  by  all  possible  collisions  is 

iM dd'  P»dD  (39) 

with  the  integration  limits  defined  by  the  range  of  size  under  consideration. 

If  the  results  obtained  in  the  previous  section  for  two  discrete  sizes  are  applied  to  the  present 
case  of  a  continuous  size  distribution,  the  quantities  in  the  above  notation  list  are  (with  the  under¬ 
standing  that  all  integrals  are  taken  from  the  minimum  diameter  to  the  maximum  diameter) 


Cn=P(D)-C 


(40) 


where  P(D)  is  the  given  probability  density  function  of  the  size  distribution, 
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■ 


Figure  6.  Shear  stress  in  a  simple  shear  flow  of  spheres  of  two  sizes. 


By  use  of  the  revised  theory  of  a  two-size  mixture,  the  comparison  between  the  theoretical  pre¬ 
diction  and  experimental  data  obtained  by  Savage  and  Sayed  (1984)  becomes 


shear  stress  (two-size  given  by  eq  37  with  n  =  3) 
shear  stress  (single-size  given  by  eq  1 3a) 


COMPLETE  SPECTRUM  ANALYSIS  FOR  SPHERES  AND  DISKS 

Based  on  the  analysis  for  a  two-size  mixture,  stresses  can  now  be  obtained  for  mixtures  with 
any  size  distribution.  A  typical  size  distribution  in  a  natural  material  is  shown  in  Figure  7a,  where 
F(D)  is  the  cumulative  volume  percentage  of  spheres  that  have  sizes  up  to  D. 

The  derivative  of  F(D)  is  denoted  by  P{D)  as  shown  in  Figure  7b.  The  volume  fraction  of 
spheres  that  have  sizes  between  D  and  D+dD  is  represented  by  P{D)dD.  The  results  obtained  pre¬ 
viously  for  a  binary  mixture  can  be  directly  transferred  to  this  complete  spectrum  analysis.  The 
two  discrete  sizes  DL  and  Ds  are  replaced  by  two  differential  size  ranges:  D  to  D+dD  and  D'  to 
D'+dD'.  The  spectra  for  both  D  and  D'  are  the  same  and  are  determined  by  the  probability  distri¬ 
bution  function  F(D)  or  the  probability  density  function  P(D). 

The  following  notations  will  be  used  in  the  succeeding  discussions: 


C 

CDdD 

GDdD 

N0dD 

P0dD 

Nbo'dD' 

^DD' 


total  volume  concentration  of  all  spheres 
volume  concentration  of  spheres  with  diameters  D  to  D+dD 
number  percentage  of  particles  with  diameters  D  to  D+dD 
number  of  spheres  in  a  unit  volume  with  diameters  D  to  D+dD 
number  of  spheres  with  diameters  D  to  D+dD  on  a  unit  surface 
frequency  of  collisions  a  D-size  sphere  receives  from  all  the  spheres 
with  diameters  D'  to  D'+dD' 

average  momentum  transfer  between  two  spheres  of  diameters 
ranging  from  D  to  D+dD  and  D'  to  D'+dD' 
average  energy  dissipation  because  of  a  collision  between  spheres 
with  diameters  D  to  D+dD  and  D'  to  D'+dD' 
fluctuation  speed  of  spheres  with  diameters  D  to  D+dD 
average  gap  between  adjacent  spheres. 


Using  the  notation  above,  we  see  that  the  rate  of  momentum  transfer  across  a  unit  surface  because 
of  collisions  between  spheres  of  diameters  D  to  D+dD  and  D'  to  D'+dD'  is 
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I  Of—  C*0.5  ^t*0.2 


Olh-  C  =  04 


0  1  0.2  0.3  0  4  0  5 

CS 


0  0.2  0.4  0.6  0.8  10 

RD 


Figure  5.  Asymptotic  solutions  for  shear  stress  in  a  simple  shear  flow  of  spheres  of  two  sizes 
(dashed  lines  =  RD  -+  0,  solid  lines  =  RD  -*  1). 


fluctuations  are  not  negligible.  However,  this  analysis  gives  a  possible  condition  under  which  Bag- 
nold’s  assumption  is  valid:  when  large  spheres  are  sheared  in  dust-like  material,  the  large  spheres 
only  move  with  the  mean  flow  and  the  resulting  stresses  have  the  same  characteristics  as  those  for¬ 
mulated  by  Bagnold  (1954).  Equation  36b  also  shows  that  as  p  increases  from  0  to  4/7,  a  reversal 
of  t2 i  from  negative  to  positive  occurs.  It  is  believed  that  if  particle  rotations  due  to  the  friction¬ 
al  impact  are  included  in  the  analysis,  such  reversal  will  be  removed  and  r22  will  stay  negative. 

In  general,  as  the  shearing  spheres  decrease  in  size,  the  stresses  decrease  too.  Equations  36a 
and  b  are  derived  for  the  case  when  RD  approaches  0  while  Cs  and  CL  are  fixed,  therefore  they 
set  up  a  lower  bound  for  the  stresses  in  a  binary  mixture.  For  a  given  set  of  Cs,  CL  and  £>L,  as 
Ds  approaches  0,  the  stresses  must  reduce  monotonically. 

Equations  13a  and  36a  are  plotted  in  Figure  5  where  the  general  behavior  of  the  shear  stress 
for  0  <  RD  <  1  is  depicted.  The  dashed  curves  are  from  eq  36a,  which  represent  the  case  as  RD 
approaches  0;  the  solid  curves  are  from  eq  13a,  which  is  applicable  as  RD  approaches  1 . 

An  analytical  formulation  of  the  stresses  for  an  arbitrary  RD  is  difficult.  However,  with  the 
two  limiting  cases  quantified  in  Figure  5,  smoothing  procedures  as  discussed  by  Churchill  and 
Usagj  (1972)  can  be  used  to  combine  the  two  curves  into  one  curve  that  describes  the  stresses  for 
the  entire  range  of  RD.  A  possible  way  to  join  the  two  limiting  cases  is  given  as: 


^RD)-T^(RD)rKx  (RD) 


P-> 

'  rii .  m  ' 


and  7jj0  are  the  stresses  given  by  eq  36a  and  b,  which  correspond  to  RD  -*  0,  and  are  the 
stresses  given  by  eq  13a  and  b,  which  correspond  to  RD  -*  1 .  The  larger  n  is,  the  closer  follows 
the  limiting  values  of  Tyo  and  Ty  ] .  Equation  37  approximates  the  stresses  for  a  two-size  mixture 
of  any  0  <RD  <1.  r21  (RD)  is  plotted  in  Figure  6  with  n  =  3.  It  was  found  that  for  n  >  3,  the 
solution  given  in  eq  37  remains  independent  of  n  up  to  the  second  decimal  place. 
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‘22 


dl  du  /dl+xl\J  c. 

'  P,T(1«)  (-0.07  *  0.04  ^ 


DS  CL  Cg  Vg  j 

^  °L  SS\ 

(31a) 

(31b) 

in  which  the  terms  involving  the  square  and  higher  powers  of  RD  are  dropped.  The  energy  dissi¬ 
pation  in  a  pair  of  like  spheres’  collision  is 


ell  ~  Ps 


2  /zV^l\2  ,1-e2  u(  1  +c)  u2(l+e)2  , 

6  Vfic/  V  2  /  1  4  7T  "  4  } 


(32) 


where  collisions  are  assumed  to  occur  with  relative  velocity  determined  strictly  by  the  mean  flow, 
and 


r  -nnDs  ,/2^2  ,  P0+e)  P2(l+e)2 

Ess-P*~vs  (4-  +  ~ —  -  — ; 4 — 


re 


-  - 

v. 

i 

I 

1 


£ 

¥ 

& 


(33) 


where  collisions  are  assumed  to  occur  at  relative  fluctuation  velocity.  Since  the  large  spheres  have 
infinite  mass  compared  with  the  small  spheres,  we  have 


^LS  Ess 


(34) 


Substituting  eq  10,  23, 27, 31a,  32, 33  and  34  into  eq  9,  we  obtain  Vg  as 


■(o^o^r-0  H 


l-g2  .  Mill)  M2(l+e)2'. 

^  T  —  ”  /i  ' 


du 

(lX  2 


(35) 


The  shear  and  normal  stresses  are  thus  obtained  after  substituting  eq  35  into  eq  31a  and  b  and 
setting Ds/Dl  equal  to  zero, 


„2/3 

lL  Co 


Tii  =  PsCDi  (^T  )*  (l+e)  (°-05  +  0.08  p)  /3  I/3 

2  7TC(C0  -  CL  ) 

-O.Q7-K).04  <2 
Tj2  0.05+0.08  p  r21  • 


(36a) 


(36b) 


An  interesting  observation  is  that  r21/r22  is  only  a  material  constant  now,  independent  of  the 
solid  concentration.  This  phenomenon  is  quite  different  from  the  result  when  RD  approaches  1 
as  shown  in  eq  13a  and  b.  Bagnold  (1954)  derived  the  stresses  assuming  that  single-size  spheres 
do  not  fluctuate  while  being  sheared  in  a  Newtonian  fluid.  His  result  also  demonstrated  that 
r21/r22  is  independent  of  the  solid  concentration.  Experiments  show  that  velocity  fluctuations 
do  exist  in  the  shearing  flow  of  a  nearly  single-size  solid-fluid  mixture  and  previous  analysis, 
which  includes  velocity  fluctuations,  shows  that  r2 , /r22  does  depend  on  solid  concentration 
(Shen  and  Ackermann  1982),  hence  Bagnold ’s  result  does  not  describe  the  flow  when  velocity 
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^Sa=  Cs/(1  -CL). 


(24) 


Let  CSo  denote  the  densest  concentration  of  the  fine  spheres, 


(25) 


After  substituting  eq  24  into  eq  25  and  solving  for  ss, 

%=^s((Cso/C’sa)‘/3 -1).  (26) 

Let  AgL  be  the  frequency  of  collisions  between  a  small  sphere  and  the  surrounding  large  spheres. 
The  total  number  of  collisions  between  large  and  small  spheres  is  equal  to  the  product  of  N^s  and 
the  number  of  large  spheres  and  also  equal  to  the  product  of  N$L  and  the  number  of  small  spheres, 
therefore, 


JVsl-M 


LS 


(27) 


When  a  small  sphere  being  considered  is  not  adjacent  to  a  large  sphere,  the  collision  frequency 
between  a  small  sphere  and  the  surrounding  small  spheres  is 

^ss  =  v$lss  ■  (28) 

If  the  small  sphere  being  considered  is  adjacent  to  a  large  sphere,  then  only  half  of  the  collisions 
of  the  small  spheres  are  with  small  ones;  the  other  half  are  with  large  ones.  However,  since  the 
number  of  small  spheres  minus  the  number  of  small  spheres  adjacent  to  large  spheres,  or 


6q 

ttZ)3 


6  Cs  6  C, 
'^1 


6  Cs  /  \ 

^3  (l  “o'  6CLU6Cs/nDl  as  Z>S/0L- 0 
St' 


the  above  result  indicates  that  as  Ds/Dt  approaches  0,  almost  all  the  small  spheres  are  surrounded 
by  only  small  spheres.  Hence,  eq  28  is  a  good  approximation  of  the  collision  frequency  between 
any  small  sphere  and  its  surrounding  small  spheres. 

As  discussed  earlier,  as  RD  approaches  0,  ^approaches  0  and  large  spheres  follow  the  mean 
flow  motion.  In  this  case  the  momentum  transfer  in  both  the  x2  and  x2  directions  in  a  collision 
of  two  large  spheres  is  a  result  of  the  relative  motion  from  the  mean  velocity  gradient  only. 

Hence, 


^LLi  =  Ps  0+<0 (0  °5  +  0.08  M) (Dh  +  sL)  d~  (29) 

which  is  the  same  as  eq  6a  in  the  previous  section,  and 

nDl  du 

^LL2  -  p>  ~r  (1+e>  (°07  +  0  04  m)  (^L  +  *L  (30) 

which  is  obtained  by  averaging  the  momentum  transfer  induced  by  the  mean  flow  in  the  x2  -direc¬ 
tion  (while  eq  6b  was  obtained  by  averaging  momentum  transfer  induced  by  the  fluctuation  in  the 
jc3  -direction).  Substituting  eq  21 , 23, 27, 28  and  29  into  eq  5,  we  obtain  the  shear  and  normal 
stresses  in  terms  of  the  unknown 
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centration  of  large  spheres  CL ,  the  average  gap  between  two  adjacent  large  spheres  can  be  com¬ 
puted  from 


Cl  _  /  A.  )3 

CLo 


where  CLo  is  the  densest  concentration  of  large  spheres.  Hence 


The  collision  frequency  between  a  large  sphere  and  its  neighboring  large  spheres,  jV£l,  is  the  ratio 
of  their  average  relative  velocity  and  their  average  gap  size.  It  is  shown  in  Appendix  A  that 


3ff/2  »/2 

l  I  £(Dl+sl) 

tt/2  0  L‘x2 _ 


cos <p  dtp  dd 


LL  3ti-/2  7t/2 

J  J  -Dl  sin0  cos 6  +  >/D2  sin2 <p  cos2  6  +  (Z)L  +  sL)2  -D\  d<j>  d6 
rr/2  0 

where  the  following  binomial  expansion  may  be  used  for  simplification, 
x/Dl  sin20cos20  +  (Z)L  +  sL)2  -D2L  =  ( DL  +  sL ) sin0  cos# 


( 1  -  sin2  0  cos2  6)(2Dlsl  +s£  )  , 

+  — r-.— -  +  -  •  (20b 

2 \/(DL+sLy  sin2 <t>  cos2 6 

When  sl/Dl  is  small,  eq  20a  reduces  to 

~  .  (21) 

LL  2  dx2  sL 

This  approximation  overestimates  N[L  if  sL/DL  is  not  very  small,  or  equivalently,  when  large 
spheres  become  sparse. 

Since  large  spheres  do  not  fluctuate,  the  collision  frequency  between  a  large  sphere  in  Figure 
4  and  the  surrounding  fine  dust-like  spheres  is 


Nls=Ps  •  al  '~2 


where  ps  is  the  number  of  fine  spheres  per  unit  area  in  the  flow  field,  Ah  is  the  surface  area  of  the 
large  sphere,  and/s  is  the  frequency  of  fluctuations  of  a  fine  sphere.  Hence 


Afa=3C. 


(Dl\2  ^ 

\  Dc  '  Sc 


where  ss  is  the  gap  between  two  adjacent  fine  spheres.  Rather  than  assuming  that  the  average  gap 
between  small  spheres  is  the  same  as  that  between  large  spheres  as  was  done  in  the  previous  section, 
a  more  realistic  model  will  be  pursued  as  follows,  which  applies  to  the  case  when  RD  approaches  0. 
Since  the  available  volume  for  the  fine  spheres  is  1  -CL ,  the  actual  concentration  of  the  fine  spheres 
that  fill  the  interstice  of  the  large  spheres  is 
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Figure  4.  Simple  shear  Jlow  of  spheres  in 
dust-like  material. 

In  the  limiting  case  as  RD  approaches  0,  it  is  no  longer  reasonable  to  assume  that  the  mean  gap 
between  small  spheres  is  the  same  as  that  between  the  large  ones.  Most  importantly,  in  this  limit 
the  equipartition  of  fluctuation  energy  defined  in  eq  12  forces  v'L  to  approach  0.  The  equipartition 
of  fluctuation  energy  still  applies  as  RD  approaches  0.  However,  the  collision  frequency  needs  to 
be  modeled  differently  in  this  case,  as  will  be  explained  below. 

As  modeled  previously,  the  collision  frequency  is  calculated  as  the  ratio  of  v  and  the  mean  gap 
size  between  spheres.  Therefore,  when  v'L  approaches  0,  the  collision  frequencies  involving  large 
spheres  approach  0.  The  stress,  which  is  now  obtained  by  collisions  between  small  spheres,  is  pro¬ 
portional  to  .  Hence  when  RD  approaches  0,  the  stresses  approach  0  (see  Figure  5a  in  Shen 
[  1985 1 ).  This  result  contradicts  reality,  however,  since  when  the  concentration  of  large  spheres  is 
high,  the  existing  mean  shear  as  shown  in  Figure  1  must  result  in  collisions  whether  particles  fluc¬ 
tuate  or  not.  In  fact,  the  collision  frequency  between  neighboring  spheres,  if  modeled  rigorously, 
should  be  a  function  of  both  the  velocity  of  fluctuation  and  the  relative  mean  velocity  ascribable 
to  the  mean  shear  flow.  It  is  when  the  fluctuation  speed  is  much  higher  than  the  relative  mean 
velocity  that  the  collision  can  be  modeled  as  a  result  of  fluctuations  only.  As  RD  approaches  0, 
v'L  approaches  0  and  the  collisions  between  large  spheres  become  a  result  of  the  relative  mean 
shear  velocity  alone. 

In  the  limiting  case  as  RD  approaches  0,  we  may  view  the  two-component  system  as  an  assem¬ 
bly  of  large  spheres  sheared  in  a  highly  energy -absorbing  material.  As  illustrated  in  Figure  4,  this 
material  consists  of  numerous  fine  dust-like  particles  that  are  fluctuating  and  spherical.  These 
fine  particles  collide  with  each  other  and  the  large  spheres.  The  large  spheres,  driven  by  the  mean 
shear  motion,  occasionally  collide  with  neighboring  large  spheres.  The  collisions  between  large 
spheres  cause  them  to  temporarily  deviate  from  the  local  mean  motion.  From  the  energy  point  of 
view,  this  temporary  deviation  of  the  large  spheres’  velocity  produces  a  fluctuation  energy  in  them. 
Because  of  the  equipartition  of  energy,  this  temporary  imbalance  of  fluctuation  energy  is  absorbed 
by  the  surrounding  fine  dust  particles.  As  RD  approaches  0  the  amount  of  energy  dissipated  per 
collision  between  the  large  and  small  spheres  remains  the  same  as  described  by  eq  1  lb  and  12.  How¬ 
ever,  as  CSICL  is  kept  constant,  and  as  RD  approaches  D,  the  number  of  the  fine  particles  surround¬ 
ing  the  large  spheres  increases  as  RD" 3 .  The  fluctuation  speed  of  the  small  spheres  and  the  colli¬ 
sion  frequencies  increase  as  RD~3h .  Hence  the  time  scale  of  absorbing  this  energy  imbalance  de¬ 
creases  as  RD 9/2  ■  Therefore,  as  RD  approaches  0,  the  large  spheres’  motion  reduces  to  the  local 
mean  velocity  long  before  the  next  collision  between  the  large  spheres  takes  place.  In  this  limiting 
case,  the  collision  frequency  between  large  and  small  spheres,  which  is  the  main  factor  of  the 
stress  generation,  is  no  longer  a  function  of  v’L,  but  rather  a  function  of  the  mean  velocity  gradient 
alone. 

Consider  randomly  distributed  large  spheres  sheared  as  shown  in  Figure  4,  with  the  interstice 
filled  with  small  spheres  with  diameter  Ds  approaching  0.  The  large  spheres  travel  strictly  with 
the  mean  flow.  The  average  collision  frequency  between  a  large  sphere  and  the  surrounding 
large  spheres  is 

K[l=VJsl  (17) 

where  FL  is  the  average  relative  velocity  between  two  adjacent  large  spheres  and  sL  is  the  average 
gap  between  two  adjacent  large  spheres.  In  a  random  assembly  of  large  spheres,  with  a  given  con- 


Figure  3.  Particle  size  distribution  of  spheres  (after  Saved  1 981 ). 


where 


RD  =  Ds/Dl  ,  RC=  CS/CL  04) 

F =  RD-1  +  /s  +  ( 1  +/,)  «C2  JUT*  +  R~R^  ( 1  +/UT*  +2/s)  (1  +/?/>' )  (1 5) 


and 


.  =RDl±R£  Cq3  - 

/s  RD3+RC  c'h  W 

We  now  compare  the  model  with  two  tests  conducted  by  Savage  and  Sayed  (1984).  Test  1 
used  nearly  uniform  plastic  beads  with  a  mean  diameter  of  1 .32  mm.  Test  2  used  a  7:3  mixture  of 
ot  the  same  material  with  mean  diameters  of  1 .65  mm  and  0.6  mm.  The  weighted  mean  of  the 
two  sizes  is  1.34  mm.  The  size  distributions  for  test  1  and  test  2  are  given  in  Figure  3  (Sayed  1981). 
The  ratio  of  the  shear  stresses  obtained  in  these  two  materials  was  measured  to  be 

shear  of  test  2/shear  of  test  1  =  1/4.5  . 

The  theoretical  result  according  to  eq  13  is 

shear  of  test  2/shear  of  test  1  =  1/10  . 

This  comparison  shows  that  the  theoretical  prediction  is  in  the  right  order  of  magnitude;  however, 
it  is  apparent  from  Figure  3  that  the  size  spectrum  in  test  1  and  the  two  individual  spectra  in  test 
2  are  not  very  narrow.  The  above  comparison  will  be  revised  after  the  analysis  for  a  continuous 
size  distribution  is  completed. 


LIMITING  CASE  OF  THE  TWO-SIZE  MIXTURE 

At  the  time  when  the  analysis  in  Part  1:  Two  Grain  Sizes  (Shen  1985)  was  carried  out,  this 
author  was  not  aware  that  the  geometry  and  kinematics  of  the  granular  assembly  as  RD  approaches 
0  is  very  different  from  the  case  when  RD  approaches  1 ,  as  will  be  discussed  in  detail  later.  An 
implicit  assumption  in  the  analysis  in  Part  I  is  that  RD  approaches  1.  Modified  constitutive  qua- 
tions  that  relax  this  assumption  are  derived  in  this  section.  These  modified  constitutive  equations 
are  obtained  by  combining  the  result  shown  in  the  previous  section,  which  applies  to  the  case 
when  RD  approaches  1 ,  and  the  result  for  the  case  as  RD  approaches  0. 
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(9) 


t21  jx^  ~  ^LS  &  LS  +  L  ^LL  +  ^SS  ^SS 

where  WPg  (P  and  Q  =  L  or  S)  is  the  frequency  of  all  collisions  between  a  P-size  and  a  0-size 
sphere  in  a  given  unit  volume  and  ^PQ  is  the  average  energy  tussipated  in  a  collision  between 
/’-size  and  0-size  spheres.  By  definition. 
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The  average  energy  dissipation  in  each  of  the  three  different  kinds  of  collisions  is  derived  as 
(Shen  1985): 


^  LL  or  SS 
£LS  = 


,2  A-€2 

~  Ps  c  VL  or  S  t  4 


rr D3,  vl'+v^  ,  l-e2 


L  ‘'S  r 
3  1  4 


Jt(l+c)  m10±£>1i 
it  '  4  J 

+  ju(l+e)  _/iLLl±f)!i 
7 r  4  J 


To  relate  v[  and  Pg  ,  an  equipartition  of  energy  is  assumed.  Then 
D[vi2  =  £>§p'2 


(Ha) 

(lib) 


(12) 


namely,  the  energy  contained  in  the  fluctuational  motion  of  the  large  and  small  spheres  is  the 
same  (the  justification  of  this  assumption  is  given  by  Shen  1985). 

After  assembling  eq  6-12  into  eq  5,  the  stresses  for  a  mixture  of  two  sizes  of  spheres  are  ob¬ 
tained  as 

=  c  d2  —  )2  (i+e)3/2(°-053+0-081  p?17  .  [i-KflOfl^r3]  Cl/3 

Ps  L^dx2  rl-€2  n(l+e)  ^(l+e)2,^  2[1  +(RC)RD~2]  (C'b -c'h) 

1  4  7T  4 


F3h  _  J _ 

[i+(/?c)  Rcr^i)2  \+RC 


(13a) 


(1+# 
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'  ‘  (0.053*0.081  ^ 


[2 VTWP-  * 

(1+/?Z?3)^ 


(13b) 
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These  three  quantities-AA/pg.,  jVpg  and  pPj-are  derived  as  (Shen  1985): 


^lli  =P,X  0+e) (0.053 +0.081  n)^  (Dl+s) 


^  ?0+ £_) 

6  7T2  ^2  \ 


AW  - t  A1 

^LLl-'Ps  6  ff: 

—  n*>s  cfu 

AMSSl  =  PS  "T  0+f) (0.053  +  0.081  ji)2L  (£>s+s) 

» rr  n®s  2(  1  +e)  /—  , 

A*SS2  ="Ps  -g - -r-  V2  vs 


7tZ)3l 


—  —  du  /As  +^*l  \  2Dg 

^sl,  =  A^ls,  =PS -^(He) (0.053 +  0.081  +*)%& 

-  *D\  4  ( 1  +e)£)|  , - - 

^SL2  =  ^LS2=-Ps  -g~  V/V^T 


SL  2 

where  v[  and  Vg  are 
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6CClZ?l 


(8a) 


6CCg  £>g 

Psi  l,l2~  n{GLDl  +  GsD3s)  ■ 


(8b) 


Since  the  fluctuations  are  caused  by  the  shearing  motion,  the  work  done  by  shear  forces  is  the 
energy  input  for  this  mode  of  motion.  On  the  other  hand,  the  collisions  produced  by  these  fluc¬ 
tuations  dissipate  energy  and  this  dissipation  is  the  energy  output  for  this  mode  of  motion.  When 
a  steady  state  is  established  in  the  simple  shear  motion  described  in  Figure  1 ,  the  energy  input  and 
output  are  balanced.  In  a  mathematical  form,  this  balance  is  presented  as 


mixture,  it  is  believed  that  although  the  results  given  in  eq  54  and  56  and  Figure  8  underestimate 
the  real  stresses  when  the  size  gradation  is  very  broad,  the  general  trend  of  stress  reduction  as  size 
gradation  broadens  is  not  going  to  be  affected. 
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APPENDIX  A:  DERIVATION  OF  COLLISION  FREQUENCY  BETWEEN  NEIGHBORING 
SPHERES  THAT  FOLLOW  THE  MEAN  SHEAR  FLOW  WITHOUT  FLUCTUATIONS 


Let  D  be  the  diameter  of  spheres,  and  s  be  the  average  gap  between  neighboring  spheres.  Figure 
A1  shows  the  location  of  a  sphere  B.  The  inner  dashed  semi-sphere  is  the  locus  of  the  centers  of 
neighboring  spheres  as  a  collision  takes  place  on  the  surface  of  B.  The  outer  dashed  semi-sphere 
is  the  locus  of  the  average  centers  of  the  neighboring  spheres  around  B. 


Figure  Al.  Geometry  of  a  collision  driven  by  mean  flow. 


The  collision  frequency  between  neighboring  spheres  is  calculated  as 

n/8  (Al) 

where  u  is  the  mean  relative  velocity  between  neighboring  spheres  and  8  is  the  me*oi  distance 
traveled  to  make  a  collision  contact.  An  example  of  a  special  8  is  given  in  Figure  Al . 

When  a  sphere  A  collides  with  sphere  B  at 


j  |  w  IT 

r  =  r  sin0  cos 6  I  +  rsin<t>  sin0  k  +  r  cos 0/  ;  -  <  6  ,  0  <  0  <  -  ,  \r\=D  (A2) 

sphere  A  must  have  traveled,  according  to  the  mean  relative  velocity,  a  certain  distance  8 /.  In  Fig¬ 
ure  Al ,  r  =  D"y  The  mean  location  of  sphere  A  before  it  started  to  move  toward  sphere  B  is,  therefore, 

R  =  (rsin0cos0  +  8)j+/-sin0sin0&+/-cos0/  ;  I/?  I  =  D  +  s.  (A3) 

Squaring  eq  A2  and  A3,  then  subtracting  one  from  the  other,  we  find  that 

8  =  -Dsin0cos0  +\/D 2  sinJ0cos20  +  (D  +  s)2  -  Di  .  (A4) 

The  velocity  of  sphere  A ,  which  has  moved  this  distance,  is 

u  =  ~(D+s)cos<t> .  (AS) 
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Substituting  eq  A4  and  AS  into  A1 ,  we  calculate  the  collision  frequency  as 


3jt/2  tt/2 

f  /  ~  (D+s)  costt>d<f>  d& 

tr/2  0  *2 


3jt/2  7t/2 

I  I  -  D  sin0  cos0  +  s/D1  sin* <p  cos2  6  +  (D  +  s)2  -  D 2  d<p  dd 
tt/2  0 


(A6) 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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